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ABSTRACT 

Advanced 3D radiative MHD simulations now reproduce many properties of the outer solar atmo¬ 
sphere. When including a domain from the convection zone into the corona, a hot chromosphere and 
corona are self-consistently maintained. Here we study two realistic models, with different simulated 
area, magnetic field strength and topology, and numerical resolution. These are compared in order to 
characterize the heating in the 3D-MHD simulations which self-consistently maintains the structure 
of the atmosphere. We analyze the heating at both large and small scales and find that heating is 
episodic and highly structured in space, but occurs along loop shaped structures, and moves along 
with the magnetic field. On large scales we find that the heating per particle is maximal near the tran¬ 
sition region and that widely distributed opposite-polarity field in the photosphere leads to a greater 
heating scale height in the corona. On smaller scales, heating is concentrated in current sheets, the 
thicknesses of which are set by the numerical resolution. Some current sheets fragment in time, this 
process occurring more readily in the higher-resolution model leading to spatially highly intermittent 
heating. The large scale heating structures are found to fade in less than about five minutes, while 
the smaller, local, heating shows time scales of the order of 2 minutes in one model and 1 minutes in 
the other, higher-resolution, model. 

Subject headings: Sun: magnetic topology - Sun: chromosphere - Sun: transition region-Sun: corona 


1. INTRODUCTION 

The last few decades have seen great progress, both 
observationally and theoretically, in understanding the 
heating of the solar chromosphere and corona. These 
advances have occure d through the Solar a nd Helio¬ 


spheric Observatory (Domingo et al. _ 1995), Transi - 


tion Reg ion and Coronal Exp lorer jHandy et ST 1999), 
Hinode ( [Kosugi et al. 2007), Solar Dy namics Obser¬ 


vator y/Atmospheric Imaging Assembly (Lem en et al. 
2012[), Interface Region Imaging Spectrograph ( [Pe Fon- 


tieu et ar]|2014 ) satellites, and ground-based observato- 


ries and instruments such as the Swedish 1-m Solar Tele¬ 
scope 


Scharmer et al. 2003 Scharmer 2006 Scharmer 


et al. 20081), the Interfe rometric Bidimensional Spectrom¬ 


eter ( iGavallinil 120061) at the Dunn Solar Telescope and 
^1 Mom oh 


GREGOR (Schmidt et al. 2012). This wealth of high 
quality data has been matched by the development of 
sophisticat ed numerical models of the outer solar atmo- 
sphere (e.g.|Gudiksen & Nordlund|2005| |Bingert & Peter 

-I y . t JK JK m mt . I ■ M i I • I I 


2011 Abbett'2007^ T'he convection zone clearly con- 


tams sufficient mechanical energy to heat the corona. 


viggo. hansteen@ast ro. uio. no 


However, the mechanisms responsible for the transport 
and release of the energy from the photosphere through 
the chromosphere to the corona remain e l usive (see, e.g. 
Walsh fc Ireland] |2003[ |Klimchuk| |2006| |Parnell fc De| 
Moortel||2012p. 

Observational constraints (|Athay fc White} 1978 |1979 


Mein fc Schmieder|1981||FQssum Carlsson|2005 ' Carls- 


son et al. 


Stem|200^ 


S ootT 

- mdi 


as well as numerical models (| Carlsson & 

icate that the power contained in acoustic 
waves is not sufficient to heat or maintain a chromosphere 
or corona as found on the Sun except perhaps in regions 
of very weak magnetic fields. Therefore, most coronal 
heating models require that the magnetic field plays an 
important role in the processes of transport, storage, and 
release of energy. 

Several models have been proposed to solve these prob¬ 
lems. These can be subdivided into two classes involving 
either wave heating, “AC” models, where the magnetic 
field is stressed at relatively high frequency compared to 
chromospheric or coronal dynamical time scales, or nano¬ 
flare “DC” models, where magnetic field stre sses are built 
up over longer time scales ( Klimchuk||2006 ). Both types 
of models face constraints set by the nature of the solar 
atmosphere. The AC models must explain the difficulties 
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of transporting wave energy through the chromosphere, 
as fast mode waves will be refracted towards regions of 
low Alfven speed. Alternately, if Alfven waves are con¬ 
sidered, the expected low dissipatio n rates of such waves 
the corona must be overcome ([Bondeson 1985 van 


jen 


Ballegooiien et ar]|2Qll[ |Asgari-lhrghi ^ van Ballegooi 


2012p. DC niodels must explain the formation of 
sutKciently large magnetic field gradients in the corona 
to give sufficient heating in the face of the intrinsically 
low coronal resistivity. For a comprehensive discussion 
about the magnetic models we refer to e.g. |Klimchuk| 
([2006|. 

in this paper we concentrate on the DC heating type 
models and consider the properties of numerical simula¬ 
tions of nano-flare heating. Thoug h several such mod¬ 
els have been presented previously dCudiksen fc Nord- 


lund 2QQ5[ [Berger fc Asgari-Targhij2UU9[|Bingert ^ Feter 


2011), these lacked the self-consistent convective motions 

included in our simulations. We believe that the models 
presented here are the first that demonstrate that the 
chromospheric and coronal magnetic field gradients built 
up through such self-consistent convection are sufficient 
to maintain coronal temperatures over an extended pe¬ 
riod. In particular, we study how the turbulent motions 
in the convection zone and photosphere drag and twist 
the magnetic field. This leads to a Poynting flux into 
the corona, with steadily increasing free magnetic energy 
stored in the field. This continuing turbulent motion in¬ 
creases the stored energy in the field until the maximum 
storage capability is reached and the magnetic gradients 
required for free energy storage no longer can be sup¬ 
ported. A quasi steady-state ensues where the energy 
injected is balanced by dissipation at small scales. 

Thus, dissipation occurs as the magnetic field attains 
sufficient gradients to generate small scale current sheets 
which, by reconnection, release the stored magnetic en- 
erg y and create a simpler magnetic topology as described 
by Galsgaard & Nordlund (1996). In general, energy 
is therefore released in the form of resistive currents, 
plasma flows, wave excitation, particle acceleration, or 
through viscosit y as plasma flows and other motions 
are th ermalized (Priest & Forbes 2000 Baumann et al. 

2M^. 


Tn the MHD type of model described here we do not 
have sufficient spatial resolution nor include enough mi¬ 
crophysics to model all these aspects of magnetic recon¬ 
nection in the corona with high fidelity. The basic as¬ 
sumption is therefore that on the scale we are model¬ 
ing the system, the microscopic description of these pro¬ 
cesses is not vital and that the large scale evolution of 
the system is insensitive to the details of how magnetic 
energy is thermalized. There is mou nting evidence that 
this assumption is not unreasonable ([Galsgaard fc Nord- 

lund|1996l|Gudiksen fc Nordlund|2(JU5| [Berger ^ Asgari 

Targhi||2009[ [Bingert & Peter||2oil[ |Pontinj|2011|). 

I'hough they do not have the spatial resolution to rule 
out that AG heating plays a significant role in coronal 
heating, the recent development of “realistic” 3D MHD 
models spanning from the p hotosphere to the corona 
(Gudiksen & Nordlund||2005 ) support the nanoflare the- 
ory. 'I'hese models show intense Joule heating events or 
the formation of current sheets throughout the chromo¬ 
sphere all the way to the corona where the Joule heating 
per particle peaks around the upper chromosphere and 


lower transition region (Hansteen et al. 2010). In sum¬ 
mary, large scale photospheric motions drive the mag¬ 
netic field gradients in the chromosphere and corona to 
small scales, thereby forcing episodic dissipation at the 
same small scales. These dissipation ev ents (nanoflares) 
were first put forward by Levine| ( 1974|) as a mechanism 
of heating the corona . The nanofi a re theory was f urther 
elaborated by [Parkerl (| 1983l |1988|) 
chnk fc Cargill "j pnOT ); Priest et al. 


Cargill (1994); Klim- 


Parker estimated an energy of about 10^^ to 10^^ J 
per event, and noted that a sufficient number of events, 
when combined over the Sun, should be able to provide 
most of the energy ne eded to explain the coronal energy 
budget (|Parker||1988|). 

Despite the success of numerical models, observat ional 
confirmation is of course critical. Several studies (|Dat- 


lowe et an|1974[ [Lin et ar]|1984[ |Dennis|[l985[ [Crosby 
et al.|199^ suggest that solar flares occur with an energy 
distribution given by the power law dN/dE ^ E~^ with 
an index a ^ 1.8, where N is the number of events and 
E is the energy of each event. However, if the nanofiares 
are to be energetica hy important to coronal heating a is 

required to be > 2 ( |Hudson|l991[). __ __ 

Several statistical s t udies (jKrucker fc Benz [1998 
Berghmans et al.||1998| IParnell ^ Jupp||20UUl jAschwa" 
Benz & Krucker 2002]) have attempted 


den et al. 2000 [ 


to compute the value of a for EUV, small-scale tran¬ 
sient brightening, or nanofiare-like heating events, across 
the solar disk, but have so far proved inconclusive. The 
models used in those studies relay upon iron lines (at 
171 A and 195 A) to d etect and diagnose the events. 
Krucker & Benz ([1998 ) developed a method to assess 
heating events by calculating the emission measure at 
all pixels for the entire time series. This model giv es a 
power-law with an in dex a that is greater than 2. |As-| 


chwanden et al. (2000) on the other hand, present a more 
restrictive method in selecting heating events; they use 
a pattern recognition code that extracts spatiotemporal 
events with significant variability. This method gives a 
power law with an index a that is less than 2. The diffi¬ 
culty in reaching an agreement on the frequency and en¬ 
ergetics of nanofiares arises from the difficulty in assess¬ 
ing when an event has happened. Without such agree¬ 
ment it is difficult to estimate whether the contribution 
from nanofiares is of the same order as the 300 W m“^ 


that is requ ired to heat the quiet corona (Withbroe & 
Noyes|[Tf7^ . 


There are also other observational puzzles that may be 
shedding light on the nature of coronal heating. Amongst 
these are the na ture of the spicule acceleratio n and heat¬ 
ing mechanism (De Pontieu et al.|2009 2011), the run of 
the differential emission measure with temperature (Pe 


ter et al.||2006 ), and the ubiquitous avera ge redshift seen 


in low er transition region emission lines (Peter & Judge 


1999). A successful model for coronal heating should be 
able to explain how these phenomena arise and predict 


strengths (e.g. 

Hansteen et al.| 

pMO Oiluri et al.|2015|). 

In this paper we expand on tJ 

[le work done by kansteen 


et al. (2010) and characterize Joule heating events in two 
models witn different initial magnetic field distributions 
and different spatial resolutions. We will discuss the 
shapes, lifetimes and sizes of the current sheets together 
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with the structure of the magnetic field lines about the 
current sheets. 

A short description of the code used and the models 
are presented in section § 2, the characterization of the 
heating events is given in § 3, and the discussion and 
conclusion follows in section § 4. 

2. NUMERICAL METHOD AND MODEL 

We wish to model coronal heating resulting from con¬ 
vective motions by consistently simulating the layers of 
the solar atmosphere from the upper convection zone to 
the corona. Solar convection is maintained by setting 
the entropy of inflowing material at the bottom bound¬ 
ary of the computational domain and solving the equa¬ 
tions of radiative transfer in four frequency bins such 
that an atmosphere with effective temperature approxi¬ 
mately equal to the so lar T^ff = 5780 K arises. This is a 
well known technique ( Nordlund||1982 Steffen & Much- 


tial magnetic field distribution. The model A described 
below is ess entially the same mod el as the ‘Bl’ model 


describ ed in Hansteen et al. (2010) and Guerreiro et al 

dMl. 


1988) that reproduces solar photospheric intensi 


ties, energetics, and dynamics to a high degree of preci¬ 
sion, seemingly only limited by the spatial resolution of 
the model. 

In order to extend the model into the chromosphere 
we include the effects of s cattering in the equations of 
radiative transfer following Skartlien (2000). 

A magnetic field must also be introduced into the 
model. This is done by setting the vertical magnetic field 
at the bottom boundary at time t = 0 s and extending 
the field into layers above by performing a potential field 
extrapolation. As the model evolves we allow the field 
to move with the fluid motions at the bottom boundary. 
In the models presented here no new flux is injected into 
the bottom boundary during the simulation run. 

The model is evolved in time by solving the radiative 
3D MHD equations using the Bifrost code. A descrip¬ 
tion of the cod e and treatment of th e MHD equations 
can be found in Gudiksen et al. (2011). In summary, the 
code includes artihcial viscous dittusivity and magnetic 
resistivity terms, which are the sources of viscous and 
magnetic heating. Non-grey, non-LTE radiative transfer 
in the photosphere, and optically and effectively thin ra¬ 
diation in the upper chromosp here are included fo l lowing 
the schemes described above (Hayek et al. 2010 Garls- 
son & Leenaarts 2012). The equation of state is given 
by a set of look up tables that given the internal energy 
and density return the pressure, temperature, opacities 
and other ra diation quantities as described in [Gndiks^ 
et al.| (2011). Thermal conduction along the magnetic 
held lines is included and is implemented using an im¬ 
plicit algorithm solved using a multi grid method. The 
boundaries are periodic in the x and y direction, while 
being non-periodic in the vertical 2 ; direction. The lower 
boundary allows flows to exit unimpeded while the en¬ 
tropy of inflowing material is set to maintain the effective 
temperature at Teff ~ 5800 K as mentioned above. The 
upper boundary is transparent. 

In both the two models (model A and model B) pre¬ 
sented in this paper we find a corona characterized by 
plasma at T > 10^ K maintained ultimately by Joule 
heating events that either heat the plasma directly, or in¬ 
duce small scale shear flows that are thermalized via vis¬ 
cous dissipation, as the initially potential magnetic field 
is stressed by convection zone and photospheric motions. 

The two models differ in the grid resolution and ini¬ 


2.1. MODEL A 

The model A has a computational box with physical di¬ 
mensions 16.6 X 8.3 X 15.5 Mm on a grid of 256 x 128 x 160 
points. The grid is uniformly spaced, about 65 km, in 
the X and y direction and non uniform in the ^ direc¬ 
tion, with Az ranging from 32 km in the photosphere to 
44 km in the corona. We have used the convention of 
setting the zero point in our height scale in the photo¬ 
sphere at the point where on average rsoo nm = 1- The 
unsigned magnetic field strength is fairly large, of order 
135 Gauss, but well mixed, with a net magnetic flux close 
to zero. Thus, this model could be typical of a strong- 
held, quiet-sun region with no held organized on larger 
than photospheric scales. 

A summary of the state of the atmosphere at t = 1000 s 
from the photosphere and up as function of height is 
shown in hgure The left panel shows the tempera¬ 
ture distribution. In the photosphere, the temperature 
varies between 5000 K and 7500 K, at ^ = 0 Mm. In 
the hrst 1000 km above the photosphere the average and 
minimum temperatures fall, while the maximum tem¬ 
perature falls to 5800 K some 350 km above the photo¬ 
sphere and increases thereafter. At 1 Mm we find that 
in regions of vigorous expansion, such as in the wake of 
shock waves, the minimum temperature falls to 2000 K. 
In order to avoid problems with the equation of state an 
artificial heating term is turned on at tem peratures lower 
than 2000 K (see Leenaarts et al. 2011). The maximum 
temperature at 1 Mm reaches 6 d 00 K while the aver¬ 
age is 4500 K. From this height and up to the transition 
region the minimum temperature remains (artificially) 
fixed near 2000 K, while the average and maximum tem¬ 
perature increase slowly. Material at transition region 
temperatures is found from 2 ; = 1.2 Mm and up to some 
3.5 Mm. Above this height material is found in the tem¬ 
perature range 5 x 10^ K to 2 x 10^ K. Note the three 
ionization bands with elliptical shape corresponding to 
the ionization of H l, He l, and He ll from bottom to 
top, respectively (see Leenaarts et al. 2011). We also 
find a large spread in the density in the upper chromo¬ 
sphere, transition region and lower corona between 1.5 
and 7 Mm as shown in the middle panel of figure 
At greater heights we find that the density is roughly 
constant; log^^o P = —15.3 g cm“^ equivalent to electron 
densities slightly less than = 10^ cm“^ in the corona. 
The right panel of figure shows the Joule heating and 
the average magnetic energy density. We note that the 
ratio of gas pressure to magnetic pressure /3 < 1 above 
some 1 Mm. The Joule heating decreases exponentially 
with height from the middle chromosphere to the corona 
with a scale height that is closely related to the magnetic 
energy (H^/2/io) scale height, roughl y ^ 1300 km in this 
model. Sini ilar findings were made by [Gudiksen fc Nord- 
lund| (|2005|) and|Hansteen et al.| (|201C)|). There is a large 
scatter m the strength of the Joule heating at any given 
height, up to 8 orders of magnitude. 

This simulation starts with a semi-random distribution 
of vertical magnetic field “poles” at the bottom bound¬ 
ary designed to give a “reasonable” unsigned magnetic 
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Fig. 1.— State of the atmosphere at t = 1000 s as function of height in model A. From left to right: temperature, density and Joule 
heating. The white dashed lines are the variables average as a function of height. The average magnetic energy is given in the rightmost 
panel as a dashed black line. 


flux and a few large scale unipolar concentrations in the 
photosphere. The initial field in the computational box 
is computed through a potential field extrapolation. The 
simulation is run for over an hour solar time, during 
which the magnetic field is stressed by convective and 
photospheric motions. The average field strength in the 
photosphere is of order 100 Gauss, starting at 150 Gauss 
at t = 0 s and falling to 100 Gauss at the end of the 
run as magnetic field is pulled down by convective mo¬ 
tions. On the other hand, in the upper chromosphere 
and in the corona the mean unsigned field strength is 
very nearly constant or even increasing; for example at a 
height of 1.5 Mm above the photosphere we find an aver¬ 
age field strength that varies in the range [55, 60] Gauss, 
while in the corona at a height of 4 Mm it is in the range 
[15,18] Gauss. 

In figurewe show the vertical field in the photosphere 
and a number of field lines extending into the corona at 
8 minutes into the run. The magnetic field configura¬ 
tion can also be recognized in the temperature structure 
of the corona where we see large loop shaped regions of 
nearly isothermal material following the magnetic field 
lines. The atmospheric structure is complex and there 
are large horizontal variations found, starting in the up¬ 
per chromosphere and transition region. 

The energy flux that ultimately heats the upper chro¬ 
mosphere, transition region and corona is mediated 
from the photosphere in the form of a Poynting flux 
(E X B) //io- The vertical component of this flux can 
be divided into two parts; one is due the work done by 
horizontal motions on the vertical field 

Bzi^yBy -{-UxBx)’) ( 1 ) 

the other due the transport of horizontal field by vertical 
motions 

^u,{Bl + Bl), ( 2 ) 

/^O ^ 

but note that this division, perhaps convenient for ob¬ 
servers, also contains contributions from the effects of 
waves and large terms of opposite sign, that will ulti¬ 
mately cancel, due to motions parallel to a tilted field. 
Energy is transported both as a result of the horizontal 
motions associated with convection and granulation and 
as horizontal field is advected in the vertical direction. 


We find that both components are important throughout 
the solar atmosphere, even when we do not introduce any 
horizontal field at the bottom boundary to model flux 
emergence. In figure we show these Poynting fluxes as 
a function of height and as a function of time at a spe¬ 
cific height 2 : = 1.3 Mm above the photosphere. While 
both magnetic fields and fluid velocities are possible to 
measure near the phot osphere, it is c lear that measuring 
the Poynting flux (see |Welsch 2015| is going to be ex¬ 
tremely difficult unless all the terms m equations and 
can be accurately estimated. In the chromosphere, tran¬ 
sition region and lower corona the two components of 
the vertical Poynting flux are highly time variable and 
largely have opposite signs: the work done by horizon¬ 
tal motions results in a largely upward directed Poynt¬ 
ing flux, while the transport of horizontal field gives a 
downward directed flux. However, the time average of 
the total vertical Poynting flux is directed upward above 
the photosphere. We find 12 kW m“^ at 300 km, at 
1 Mm this has fallen by almost an order of magnitude, 
to 1.5 kW m“^, and at 3.5 Mm is further reduced to 
500 W m“^. Only some 100 W m“^ remains at 10 Mm. 
The right panel of figure shows that after an initial 
transient, the total Poynting flux varies on a time scale 
of roughly 5 minutes and is dominated by chromospheric 
and p-mode like oscillations in the atmosphere. During 
the first 15 minutes the Poynting flux is large, as stresses 
are built up in the magnetic field, while the Joule dissipa¬ 
tion slowly rises as the field gradients approach a critical 
value at which dissipation is balanced by the buildup of 
gradients. After this time the total Joule dissipation is 
remarkably constant and well balanced by the total ra¬ 
diative losses. We shall see that on smaller scales the 
Joule dissipation is much more variable, but on large 
scales the system settles fairly rapidly to a quasi steady 
heating regime. We note that there is systematic sub¬ 
duct ion of magnetic field during this run. Thus, on time 
scales longer than the runs discussed here chromospheric 
and coronal heating will gradually be reduced unless the 
magnetic field is replaced, either through a local surface 
dynamo or through the emergence of new flux from the 
deeper convection zone. 


2.2. MODEL B 
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Fig. 2. — Magnetic field configuration shown for loops with apex at different heights (top panels: side view-left, top view-right). Temper¬ 
ature stratification (bottom panels: side view-left, top view-right). The field lines have been given different colors according to their apex 
height above the photosphere. The temperature color scale is green and turquoise for the upper chromosphere and low transition region 
temperature, blue for the upper transition region temperature, purple for the low corona and red for the hottest coronal regions. The Bz 
magnetic field component, (black negative polarity, white positive polarity) is shown in the photosphere at 2 ; = 0 Mm where rsoo = 1- 
These representations are for Model A. 


Model B has a computational box with physical dimen¬ 
sions of 24 X 24 X 16.8 Mm in a grid of 768 x 768 x 768 
points. The grid extends approximately 2.4 Mm below 
^500 = 1 and 14.4 Mm above. Similar to Model A the 
grid is uniformly spaced in the x and y direction, about 
31 km and non-uniform in the 2 : direction with a mini¬ 
mum spacing of 12 km in the photosphere up to 82 km in 
the upper corona. The experiment has been run at this 
resolution for about 30 minutes solar time, after hav¬ 
ing been run at lower resolution for an extended period 
before that. It then takes of order 10 minutes for the 
simulation to adjust to the new resolution. We defined 
t = 0 s as the instant in which the simulation was in¬ 


terpolated for the present resolution of model B. This 
model is constructed to have two more or less unipolar 
magnetic regions separated by some 8 Mm. The aver¬ 
age unsigned field strength is 50 Gauss and the net field 
strength is close to zero. This model is designed to be 
typical of a quiet sun network region. 

The distribution of the density, temperature, and heat¬ 
ing as a function of height in this model is roughly similar 
to Model A, as can be seen in figure The temper¬ 
atures in the left panel show a minimum temperature 
of 2 X 10^ K and a maximum of 1.6 x 10^ K. The upper 
chromosphere and transition region are even more spread 
out in height than in Model A and we find relatively low 
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Fig. 3. — Left panel: The Poynting flux as a function of height every 200 s from t = 1000 s to t = 3600 s (colored lines from purple 
through blue, turquoise, yellow to red) and time average (thick black line) for the low resolution model. Right panel: Poynting flux at 
2; = 1.3 Mm divided by the distance of this height to the top of the simulation, along with the total spatially averaged Joule dissipation 
rjj^ (dotted line) and total spatially averaged radiative losses (dot-dashed line) above this height. In both panels the dashed lines show the 
vertical Poynting flux associated with horizontal motions, while the solid lines show the total vertical Poynting flux. These plots use data 
from model A. 


temperatures and high densities to cover the height range 
between 2 Mm and 6.5 Mm. The average Joule heating 
with height also here in Model B has a scale height that 
is similar to the scale height of the magnetic energy den¬ 
sity, and is of the order of ^ 1600 km. The Joule heating 
above 2 Mm has a scatter ranging roughly 8 orders of 
magnitude. 

The bipolar structure of the magnetic field is evident in 
the photosphere as two regions of concentrated magnetic 
field of opposite polarity, as shown in figure The Bz 
component of the magnetic field in the photosphere at 
z = 0 Mm is shown from the side (left panels) and from 
the top (right panels). Magnetic field lines are drawn to 
give an impression of the topology of the stronger parts 
of the field and are seen to stretch between regions of op¬ 
posite polarity reaching various heights. This particular 
figure shows the situation at 23.8 minutes after the start 
of Model B run. 

Aspects of the temperature structure are shown in the 
bottom panels. Upper chromospheric temperatures are 
largely confined to heights around 2 Mm, but one also 
sees isolated cool loops stretching much higher. The 
same is true for the lower transition region temperatures. 
Coronal temperatures are found everywhere above 2 Mm, 
tending to be found almost universally at the greatest 
heights, but we find that the hottest isolated loops of¬ 
ten are low lying and stretched between regions of the 
strongest field concentrations. 

2.3. GLOBAL DISSIPATION PROPERTIES 

As shown in figure and figure in the upper chromo¬ 
sphere and above, the Joule dissipation in these models 
has a scale height that is roughly the same as that of the 
magnetic energy, B‘^ /2fiQ. We find a scale height of some 
650 km above 2 ; = 500 km in Model A and z = 900 km 


in Model B. Below this height, in the photosphere and 
lower to mid chromosphere, the heating scale height is 
much smaller, only 100 km in Model A and 70 km in 
Model B. 

In the mid to upper chromosphere and above, the 
heating scale height is larger than the 200 km pressure 
and/or density scale height in the chromosphere, but is 
much smaller than the coronal scale height of 50 000 km. 


Therefore, as shown by 


Hansteen et al. 


(2010), the Joule 


dissipation per unit mass (i^j^/p), or equivalently per 
particle, show a maximum in the transition region and 
low corona: in the chromosphere the number of parti¬ 
cles decreases very rapidly and the heating per parti¬ 
cle grows exponentially with height, eventually raising 
temperatures to coronal values. At this point the pres¬ 
sure scale height is much larger than the dissipation scale 
height and the heating per particle decreases, also expo¬ 
nentially, but with the longer scale height of the magnetic 
energy density. This is shown in figure along with the 
average temperature, as a function of height, showing the 
situation at three different times for both Model A and 
Model B. 

The maximum average Joule dissipation per unit mass 
is about one order of magnitude higher for Model A than 
for Model B, presumably due to the stronger average field 
in Model A (a factor x3), but perhaps also due to the dif¬ 
ferent field topologies of the models. The decrease of the 
heating per particle with height is less steep for Model 
B than for Model A. The magnetic field scale height for 
Model A is smaller than the magnetic field scale height 
for Model B above the height of maximum Joule dis¬ 
sipation, which is placed above 2 Mm in both models. 
Note also that the average temperature rise in Model 
B is less steep than in Model A, a consequence of the 
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Fig. 4. — Same as Figurebut for Model B. 


greater amount of cool material at great height in Model 
B simulation. 

As should be obvious from figures and the Joule 
heating is not all found in events of a certain magnitude, 
events are rather spread in energy over several decades. 
It is difficult to isolate events in space and time, since 
events move and overlap, but we can make a distribu¬ 
tion function of heating magnitudes at a given height. 
In figure we plot histograms of the magnitude of Joule 
heating in individual mesh points in a 200 km high zone 
centered on z = 2 Mm at various times for Model A 
and B using a bin size of 0.01 dex in energy. Initially, 
both models show very few grid cells with large heat¬ 
ing but within a few hundred seconds a more or less 
constant distribution is established with a given slope. 
Since there are a finite number of computational cells 
the distribution is necessarily truncated at low energies 
and a peak appears; at 10“^ W m“^ for the Model 
A, and at 3 x 10“^ W m“^ for model B. However, the 
slopes should be taken with a grain of salt since the cho¬ 
sen bin size is an important parameter in the resulting 
slope: increasing the bin size would make a steeper dis¬ 
tribution while a smaller bin size would result in flatter 
distributions. There is no physical argument to choose 
one bin size over another and the index of the power 
law is not a reliable measure of the role of low energy 
versus high energy events in the corona. On the other 
hand figure does show that these models establish a 
steady state heating within a thousand seconds with a 
given spread in energetic and less energetic heating. 

3. JOULE DISSIPATION CHARACTERIZATION 
3.1. HEATING EVENTS 


On the largest, global, scale we find that the Joule 
heating in the model is (very) roughly proportional to 
the magnetic energy density. With the magnetic field 
topologies considered in this paper this implies that the 
heating is vertically stratified with a scale height of some 
thousands of kilometers, and that the heating per parti¬ 
cle therefore is concentrated in the upper chromosphere, 
transition region, and lower corona as described in the 
previous section. Let us now consider the heating on 
smaller and horizontal scales. As the (originally poten¬ 
tial) magnetic field is stressed by convective and granular 
motions, gradients build up. These gradients rapidl y lead 
to the formation of current sheets as described by Gals-| 


gaard & Nordlund (1996) along which dissipation and 
hence heating occur. 

The spatial distribution of the current sheets is quite 
complex, as can be seen in figure where the current 
density per particle in both models is shown. Note that 
though quite intermittent in both models, the current 
sheets generally follow the global structure of the mag¬ 
netic field as field lines loop from one photospheric polar¬ 
ity to the other. The complexity of the photospheric field 
in Model A is somewhat greater than Model B: positive 
and negative polarities are quite intermixed in Model A 
and there is also a larger net signed magnetic field so that 
a significant number of field lines escape the top bound¬ 
ary of Model A. The highly intermixed polarities imply 
a smaller magnetic energy density scale height, and in¬ 
deed we find that the strong currents sheets in this model 
are concentrated at low heights outlining loop shaped 
structures. In addition, we find regions of strong cur¬ 
rents that extend significantly higher than the transition 
region along the “open” field lines that pierce the top 
boundary. Model B has a smaller unsigned flux and the 
opposite magnetic field polarities in the photosphere are 
better separated. Thus, we find a larger magnetic field 
energy density scale height and fewer “open” field lines: 
regions of strong current sheets extend higher in Model 
B, but we do not find “fingers” of dissipation extending 
towards the upper boundary. It also appears that strong 
current sheets are more spatially intermittent in Model 
B. 


While current sheets stretch out, forming arc or finger 
shaped structures along the field, they rapidly collapse to 
dimensions of only a few grid zones perpendicular to the 
field. This latter scaling is ensured by the functional form 
of the resistivity which is constructed proportional to the 
grid size As, resulting in the grid Reynolds numbers of 
order slightly greater than one that are required in order 
to resolve the current sheets on the chosen grid. In the 
two simulations described here this gives a current sheet 
thickness of some 120 km for Model A and some 80 km for 
Model B. Of course, on the Sun, the current sheet thick¬ 
ness in the corona is expected to be much, much smaller, 
perhaps on the order of centimeters or less depending 
on the (unknown) physical process that ultimately halts 
current sheet collapse. Ne vertheless, as has been pointed 


out by numerous authors dGalsgaard & Nordlund] 

1996 

Nordlund & Galsgaard|1997[ 

IGudiksen Nordlund 
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Fig. 5. — Magnetic field configuration shown for loops with apex at different heights (top panels: side view-left, top view-right). Tempera¬ 
ture stratification (bottom panels: side view-left, top view-right). The temperature color scale is green and turquoise for the chromosphere, 
blue for the upper transition region, purple for the low coronal temperature and red for the hottest coronal regions. These representations 
use data from Model B. 


Bingert & Peter 2011), the total dissipation in a given 
current sheet is independent of the resolution As. On 
the other hand, while the total dissipation remains the 
same, the scale at which the dissipation takes place will 
change with resolution, and this change may have an im¬ 
pact on the local observable consequences of the energy 
dissipation. 

Following Baumann et al. (2013), the current density 
in the current sheet is given by 


jc = V X B ^ 


AB 

As 


sin (l)Bcs 
As 


(3) 


where (p is the angle characterizing the difference of di¬ 
rection of field lines on either side of the current sheet. 


set by photospheric motions, and Bq^ is a typical field 
strength just outside the current sheet. For the magnetic 
topologies considered in this paper and in the absence of 
flux emergence we typically find small angles p in the 
current sheets that dissipate the stresses passed to the 
outer atmosphere from photospheric motions. A typical 
example of this is shown in figure which shows details 
of a dissipating current sheet at t = 500 s in Model A. In 
the upper two panels we show the current sheet as seen 
from the side and from above — note that the current 
sheet shows some internal structure as it is beginning to 
fragment. Selected field lines passing above and below 
the current sheet are drawn showing the (small) angle p 
that induces the current. We note that the small angles 
found here for the onset of heating episodes (and later. 
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Fig. 6. — Average Joule heating per unit mass [solid line], and average temperature [dashed line] versus height ai t = 1500 s (blue), 

t = 2500 s (green) and t = 3600 s (red) for Model A (left) and at t = 150 s (blue), t = 950 s (green) and t = 1750 s (red) for Model B 
(right). 



Fig. 7.— Histograms of the Joule heating in a region roughly at 2 Mm height for Model A (left) and Model B (right) using a bin size of 
0.01 in the log. The dotted lines have slope equal to 1.75 for Model A and 1.3 for Model B and are adjusted by eye. The colors represent 
the times of the different snapshots: from purple in the beginning of the run to red at the end. In Model A we plot the distribution at 
t = 10 s, 100 s, 150 s, 200 s, 300 s, 400 s, 500 s, 1000 s, and thereafter every 500 s until 3500 s, and at t = 10 s, 50 s, 100 s, 150 s, 200 s, 
600 s, and thereafter every 200 s until 2000 s for Model B. 


for Model B) are not consistent w ith the “secondly in- 
stability” (SI) advocated by {e.g. Dahlburg et al. 2005) 
which appear to require substantially larger shear. 

The time evolution of this (typical) current sheet is 
shown in figure pT] where the current density squared per 
particle is shown over 900 s at intervals of 100 s. Ini¬ 
tially, there is little current density and the direction of 
field lines is nearly the same above and below the forming 
current sheet. After a few hundred seconds the differen¬ 
tial angle has increased markedly as has the current den¬ 
sity and a region of strong dissipation has formed. The 
dissipation region does not stay in the same place, nor is 
it necessarily tied to the same field lines, but moves up¬ 


ward in the atmosphere as well as horizontally (upwards 
in the series of panels shown). The upward displacement 
is of the order of 1 Mm. The horizontal displacement 
is most clearly visible from t = 500 s onwards as the 
strongest part of the current sheet moves horizontally 
^1.6 Mm with a velocity of ^ 4 km s“^ (this motion 
may be the equivalen t of that reported observationally by 
Schrijver et al. || 1999 ). Note that the differential angle in 
the magnetic field decreases as the current sheet moves 
away from the region marked by field lines. Towards 
the end of the time series the current sheet, now located 
higher in the atmosphere and above the region marked 
by magnetic field lines in the figure, has dissipated. At 
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Fig. 8.— 
view-right]. 


Current density squared per particle (j^/p) (arbitrary units), for Model A (top) and Model B (bottom) [side view-left, top 


the same time a new current sheet seems to be forming 
in the location originally stressed; the magnetic field is 
rooted in the photosphere and convection zone below, in 
which the typical time scale of plasma motions stressing 
the field are much longer than those found in the chro¬ 
mosphere and corona, allowing the repeated rebuilding 
of stress in a given upper atmospheric location. 

At certain locations, such as the one described above. 
Model A shows recurrent behavior, or periodicities, in 
the dissipation. This is shown in figure [m We select 
a fixed volume of 1 x 1 x 1 Mm^ around the top of the 
current sheet studied and integrate the total heating in 
this volume for the entire time series. Dissipation is ini¬ 
tially at a very low level as photospheric motions have 


not yet had time to build up sufficient stresses in the 
chromospheric and coronal field. After some 10 minutes 
we find a period of strong (and repeated) dissipation; this 
is the same series of current sheet formation and dissipa¬ 
tion that was described above in figures and jT After 
this series of large events we find smaller maxima after 
30 and 45 minutes, before a very large dissipation event 
after 53 minutes. Each individual dissipation event lasts 
some few hundred seconds. Both complex and simpler 
temporal structures of the maxima are found, with ei¬ 
ther a single maximum or repeated maxima separated in 
time by some hundred seconds. 

In Model A the magnetic topology is such that we find 
that a number of field lines are open and pierce the up- 
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Fig. 9.— Topology of the magnetic field in a region of strong Joule heating per particle at t = 500 s for Model A; [left - side view, 
right-top view] and current density squared per particle ((j^/p)), the color table is the same as in Figurej^ The magnetic field loops have 
three colors to distinguish the loops above (blue), very close to (red), and below (green) the dissipation region. The right panels have a 
4.4 X 3.9 Mm field of view. 


per boundary. In addition, there are also several field 
lines that are “quasi-open” in the sense that they pass 
through the upper corona several times as a result of 
the horizontal periodic boundary condition forming very 
long loops before returning to the lower atmosphere and 
photosphere. 

As mentioned above we also find significant heat¬ 
ing at these locations that could be similar to t he fan 
shaped quasi-sepa ratrix layers (QSL) regions (see [Priest 


fc Demonlin 1995|) seen i n TRA CE observations as dis¬ 


cussed byjbcnrijver et al.|(|2010|). In figure fT^ the current 
density per particle is shown tor such a region, at time 
t = 500 s, showing a “finger” of large heating extending 
to some 7 Mm above the photosphere, much higher than 
for the closed loops found in regions where the field is 
more horizontally oriented. When visualizing the mag¬ 
netic field lines in the vicinity of the heating event we see 
that the field appears twisted, and in addition, that at 
the greatest heights the field spreads out in a fan shaped 


structure. Similar to what we found for the closed loop, 
the heated region moves of order 1 Mm during a 1000 s 
period and the heating remains high as long as the field 
lines are twisted; when the twist dissipates the heating is 
reduced regardless of whether or not we see a fan shaped 
structure in the field. 

The heating evolution of this region as a function of 
time in a 1 X 1 X 1 Mm^ volume centered on the twist 
seen in figure is shown in figure 

Heating events in Model B show the same general ten¬ 
dency as found in Model A: heating events in closed 
loops are associated with small angle changes in the 
magnetic field direction that last of order 100 s before 
they dissipate as the field becomes more nearly paral¬ 
lel. The location of greatest heating moves with the 
general plasma flow as field lines are advected by the 
photospheric flow. In most, if not all cases, we do not 
see any large changes in the topology of the field or ob¬ 
vious bi-directional jets associated with the reconnection 
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Fig. 10.— Magnetic field evolution at the Joule dissipation location shown in Figure!^ at t = 100 s (top left), t = 200 s, t = 300 s (top 
right), t = 400 s (middle left), t = 500 s, t = 600 s (middle right), t = 700 s (bottom lefej, t = 800 s and t = 900 s (bottom right) from a 
top view. The panels display the evolution of the tilt in the magnetic loops and the horizontal displacement of the dissipation. The arrow 
points at the current sheet of interest. The panels have a 5.5 x 5.0 Mm field of view. These representations use data from Model A. 



Fig. 11.— Heating evolution in two boxes, 1x1x1 Mm^ (solid) 
and ~ 0.2 X ~ 0.2x ~ 0.2 Mm^ (dashed) at the top of the loop 
shown in Figure Data from Model A. 


that dissipates the stresses built up by photospheric mo¬ 
tions. The field in long, or open, loops that are strongly 
heated shows evidence of twist and often spreads out in 
fan shaped structures. 

3.2. LOCAL DISSIPATION EVENTS 

Let us now consider heating events at specific locations 
in the region of most effective heating per particle, i.e., 
in the vicinity of upper chromosphere, transition region, 
and lower corona. As already alluded to in figures pT] and 
[^the heating in a given location will be episodic, but 
sometimes recurring. We have searched for sites of great¬ 
est heating per particle, and in figures \T ^ and [T5] we show 
three distinct locations where vigorous dissipation occurs 
in Model A and B, respectively. We selected these loca¬ 
tions by averaging the Joule dissipation per particle over 
the entire time series and choosing the locations (cells) 
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Fig. 12. — Magnetic field topology surrounding a dissipative “fin¬ 
ger” extending well into the corona along a long “open” loop in the 
low resolution model at t = 500 s. The current density squared p^ 
particle ((j^/p)) is shown using the same color table as in Figure!^ 
Also shown are the approximate locations of the two fixed box^, 
size 1x1x1 Mm solid line, size 0.2 x 0.2 x 0.2 Mm da shed line) 
used to follow the heating evolution in area in Figure Note the 
twist lower down, and the fan-like structure towards the top. The 
right panel field of view is 3.0 X 2.7 Mm. Data from Model A. 


with the highest average, excluding locations neighboring 
other strong heating cells. 

The dissipation per particle in fixed volumes of roughly 
195 X 195 X ^ 200 km centered on the grid cell with the 
highest average value are then found. The same method¬ 
ology was used for both Model A and B. Even though the 
volumes integrated are roughly the same in both models 
the results are significantly different as discussed below. 

In model A, the three strongest events are all found at 
heights between 2 ; = 1.5 Mm and 3 Mm, on average at 
2 Mm above the photosphere. While all locations show 



0 10 20 30 40 50 60 

time [min.] 

Fig. 13. — Heating evolution in a 1 x 1 x 1 Mm^ volume (solid) 
and ~ 0.2X ~ 0.2x ^ 0.2 Mm^ (dashed) at the center of the twist 
as shown in Figure Data from Model A. 


large variations in the heating rate, all have significant 
heating during most of the model run. As in the exam¬ 
ples of the closed loop discussed in the previous section 
the heating rate shows peaks lasting some 100 — 200 s 
that recur on average 9 times during the hour, or every 
6 — 7 minutes, though with varying peak intensity. The 
lifetime as measured in this way is clearly a combina¬ 
tion of the co-moving lifetime of the event and the time 
it takes the dissipation region to cross the measurement 
volume. For reference, the average plasma velocity in 
this region of the atmosphere, 2 Mm above the photo¬ 
sphere, is of order 10 km s“^ so a 200 km wide structure 
should take roughly 30 s to pass through a volume with 
the dimensions chosen here. 

Figures [T^ and also display the velocity along the 
magnetic ndd, and in the lower panels, the temperature 
and the density. The velocity along the magnetic field 
shows peaks that seem to be correlated with the heat¬ 
ing events, presumably caused by the rapid expansion 
caused by a drastic heat input, but possibly also due the 
sudden pairing of two flux tubes with significantly dif¬ 
ferent thermodynamic properties. The temperature and 
density also seem to be correlated with the vigorous heat¬ 
ing events, however the correlation is not nearly as clear 
as for the velocity, and while the temperature generally 
increases with increased heating, the density can either 
increase or decrease in response to an increased heating 
rate. 

While large heating per particle in Model A is concen¬ 
trated near 2 Mm, in Model B the average height of the 
largest events is found to be nearly 3.5 Mm, spanning a 
range from 2.9 Mm to 3.9 Mm. Heating events in Model 
B also appear much more sporadic than in Model A and 
the peaks are narrower in time, as can be seen in fig¬ 
ure where the heating is displayed in a logarithmic 
scale. On the other hand the peak values are compara¬ 
ble, 50 MW kg“^ in both cases. As noted earlier, the 
heating in Model B is quite intermittent and therefore 
there is no evident recurrence of heating at the locations 
of the volumes chosen. This large intermittency could 
be due to increased fragmentation of current sheets al- 
lowed by the better spatial resolution of Model B (see 
Galsgaard & Nordlun(i||1996|). 

T'he velocity, temperature and density behavior is sim¬ 
ilar to that shown in Model A. The maximum absolute 
values of the velocity along the line seems to show a good 
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Fig. 14. — [Top] Joule heating per unit mass (solid line) and velocity along the magnetic field (long dashed line), [bottom] temperature 
(solid line) and particle density (long dashed line) as a function of time for the entire run of Model A at three locations with high dissipation 
in the low corona/transition region. 
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Fig. 15. — Same as Figurebut for Model B. 
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correlation with heating events in these locations. 

A somewhat better estimate of heating event lifetimes 
in a given location can be found by calculating the auto¬ 
covariance of the dissipation time series at that location. 
We can then use the derived auto-covariance against the 
lag to determine the duration of the dissipative events 
in each cell. Here we define the lifetime of a dissipative 
event at a given location to be equal to the width of the 
auto-covariance curve at half maximum. We considered 
a lag of up to 400 s which is longer than the estimated 
lifetime of the large scale dissipative structures. This 
choice of lags proved to be reasonable since the number 
of localized heating events that last longer than 400 s is 
negligible in our models. 

Figures and show maps of the auto-covariance 
time scales at four different heights in the atmosphere 
for Model A and B, respectively: in the low chromo¬ 
sphere (490 km), in the middle chromosphere (1 Mm), 
in the transition region/lower corona (3 Mm), and in the 
corona (9.2 Mm). These time scales are computed us¬ 
ing a cell size of 65 x 65 x ^ 40 km^ for Model A and 
31 X 31 X ^50 km^ for Model B. Also shown are the aver¬ 
age heating rates per particle and the unsigned magnetic 
field strength at t = 2000 s for the same heights. 

The lower chromospheric auto-covariance time scale 
map is very highly structured and does not bear any 
clear relation to the map of average Joule heating per 
unit mass in the same region. The time scale for heat¬ 
ing events at this height is of order a minute or longer 
(see below). Regions of long-lived heating form roughly 
elliptical or circular structures with dimensions of order 
1 Mm or smaller. The map of the time scales in the mid¬ 
dle chromosphere is very similar to that found 500 km 
below, though at this height the elliptical long duration 
structures are perhaps somewhat more eccentric. Again, 
there is no clear correlation between the time scale map 
and the average heating per unit mass map. Perhaps 
this is due to the dominance of acoustic shocks at these 
high-to-medium plasma-/3 heights, where wave motions 
push and deform fields, increasing the Joule heating rate 
as oppositely directed fields are brought together at the 
whim of plasma dynamics. 

In the low plasma-/3 transition region, 3 Mm above the 
photosphere, we see a marked change in the structure of 
the auto-covarience time scale map: the general impres¬ 
sion is that structures are more ordered by the magnetic 
field and that long duration regions are larger. These 
regions are clearly more oblong at this height and now 
bear some relation to the structure defined by the aver¬ 
age heating rate. In the corona, at 9.2 Mm, long duration 
structures mainly form linear shapes with lengths of up 
to at least 5 Mm, though there is also a large nearly 
circular structure visible. The relation to the average 
heating map which shows similar linear structures in the 
same locations is quite evident. 

The time scale maps of Model B, figure pT| are quite 
similar to those found for Model A, but in general we find 
that Model B is characterized by shorter time scales and, 
of course, more highly structured shapes. This could be 
due to the higher resolution in the Model B simulation 
that allows the heating events to fragment into intermit¬ 
tent smaller scale structures. At this resolution the cor¬ 
relation between the average heating map and the time 


scale map seems better, extending even down to the mid¬ 
dle chromosphere. In the transition region and corona 
the relation is quite obvious: regions of large heating are 
also regions of longer lifetimes. Again, we find circu¬ 
lar or slightly elliptical shapes in the chromosphere, and 
much more linear structures in the transition region and 
corona. 

As discussed above it should be clear that the timescale 
derived will be a combination of the timescale for a given 
structure to cross our chosen volume and the intrinsic 
co-moving lifetime of the heating event itself (very short 
timescales can occur if there is very little variation in the 
heating rate — in that case the auto-covariance timescale 
is determined by the small random fluctuations), possi¬ 
bly modified by the occurrence of two or more heating 
regions in the same volume. We can shed some light on 
this issue by varying the size of the grid cells utilized 
and repeating the analysis described above. Since we are 
computing the auto-covariance for each cell in the sim¬ 
ulated atmosphere there are cells where the variations 
of the heating are very small; i.e. “background cells”. 
This implies a very slow decrease of the auto covariance 
against the lag curve. Since the width at half height 
of this curve is the parameter that gives the lifetime of 
the events, these background cells will contribute with 
an overestimate of the heating events lifetime. The av¬ 
erage values found are therefore an upper limit for event 
lifetimes. 

We have rebinned the computational box with cells of 
different dimensions, keeping the same physical dimen¬ 
sion for both models. Figure displays the average 
“width at half height” of the co-variance function, i. e the 
dominant lifetime, versus height for the two models, for 
cells of different physical dimensions. For both models 
the lifetime increases with physical dimension, presum¬ 
ably as the volume considered grows to encompass the 
region containing the entire lifetime of any heating event 
and/or several heating events. It is clear that in Model 
A the average time versus height is greater for all cell 
dimensions and at all heights, though the difference is 
smaller in the vicinity of the photosphere. 

Based on figure [Ts] we find that dissipative events have 
lifetimes between 2 and 7 minutes in Model A and be¬ 
tween 1 and 5 minutes in Model B. If there is only one 
event per considered box size then these lifetimes should 
give upper limits to the total co-moving lifetime of a 
single event. This estimate assumes that, using a typical 
speed of < 10 km s“^, dissipative events cannot move fur¬ 
ther than 4 Mm (4200 km) in 7 minutes. Thus, the time 
scale found when considering cells of 4 x 4 x 4 Mm^ should 
be the greatest possible timescale for a single event. This 
gives upper limits of co-moving dissipation lifetimes of 
7 minutes for Model A and 4 minutes for Model B, in 
the transition region and corona. Examining even larger 
boxes would imply that one is looking at several events 
and the timescale will increase accordingly, asymptoti¬ 
cally tending towards infinity as the number of events 
grows, this is most probably the case for the largest boxes 
of 4 X 4 X 4 Mm^. 

4. DISCUSSION AND CONCLUSION 

We present two models with fully self consistent coro¬ 
nal heating, and analyze the characteristics of this heat¬ 
ing in the transition region and corona. The velocity 
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Fig. 16.— Magnetic field [G] at t = 2000 s (left), average Joule dissipation over the time series [kW kg (middle) and logarithm of the 
auto-covariance time scale [s] (right) at 2 ; = 0.490, 2 ; = 1.002, 2 ; = 2.992 and 2 ; = 9.244 Mm (from the top to the bottom) for Model A. 


pattern in the photosphere — granulation and associ¬ 
ated flows — is driven by self-consistent convection, mod¬ 
eled to several megameters below the photosphere, rather 
than being given by a s tatisti c al description as done by 
Gudiksen & Nordlundj (|2005|); |Bingert & Peter| (|2011|). 
in addition the models discussed here have signiticant 
improvements in spatial resolution compared with these 
previous models. As time passes, magnetic fields are 
braided by photospheric and convective motions and gra¬ 
dients build up in the chromospheric and coronal fields in 
the form of current sheets. We find significant Joule heat¬ 
ing in the vicinity of the current sheets, and this heating 
comes to dominate the energetics in the low plasma-/d re¬ 
gions of the atmosphere, i.e. at heights greater than some 
1 Mm above the photosphere up into the high corona. 

The general pattern of heating found in the two models 
presented here is very similar, but Model B has signifi¬ 
cantly greater heating in the photosphere and somewhat 
greater in the upper atmosphere, even though the av¬ 
erage unsigned magnetic field is roughly the same. In 
both models the heating we find is roughly proportional 
to the magnetic field energy density, but the topology 
of the field likely also plays an as yet not wholly deter¬ 
mined role. The distribution of the magnetic field in the 
photosphere determines the magnetic field scale height, 
with closely spaced magnetic polarities of opposite sign 
giving small scale heights and widely separated polarities 
of opposite sign giving larger scale heights. This in turn 
determines the location of the maximum heating per par¬ 
ticle, given by the ratio of the magnetic and plasma scale 
heights in the chromosphere. The fall off of heating per 


particle with height in the corona will also be smaller 
the greater the magnetic scale height; based on this ar¬ 
gument we expect unipolar plage regions to still peak in 
the upper chromosphere and transition region, but also 
have high heating rates extending to large heights. The 
results presented here are therefore probably most repre¬ 
sentative of the coronae above the quiet sun or enhanced 
network. For these “normal field configuration” regions, 
the simulations show that the Joule heating events are 
most effective in the vicinity of the upper chromosphere, 
transition region, and lower corona. 

On average we find that the corona takes a long time 
to heat and fill with mass, and the question could be 
raised as to whether a coronal loop ever reaches an equi¬ 
librium, i.e. that the coronal mass is consistent with the 


space varying 

, impulsive heating 

: (see e.g. 

Bradshaw & 

Gargill 2010 

Warren et al. 

2011| 

). The total energy in- 


scales simulated here. On the other hand, the heating is 
spatially intermittent and temporally episodic. Likewise, 
the Poynting flux injected into the corona varies strongly 
with time, but on average gives a consistent, and as com¬ 
pared to estimates based on observations, reasonable en¬ 
ergy flux into the corona. The highly variable nature 
of the Poynting flux and the necessity of measuring all 
the vector components of the velocity and magnetic field 
indicates that it may be difficult to measure the heat 
flux injected into the uppe r atmosphere based on direct 
photospheric observations (|Welsc h|2Q 15[). 


Spatially, the heating is concentrated in current sheets 
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Fig. 17.— Same as figure [T^ but for Model B. 
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the thicknesses of which are set by the numerical reso¬ 
lution, that stretch along loops. Some current sheets 
fragment with time, spawning smaller current sheets as 
they are driven to smaller scales by photospheric forc¬ 
ing motions: i.e. small heating events are aligned in 
long strands, oriented along the magnetic field lines and 
associated with small changes in the angle of the mag¬ 
netic field over the current sheet where heating is con¬ 
centrated. These strands therefore show loop like struc¬ 
ture with lifetimes of the order of 100 s. This tendency 
towards fragmentation and intermittent heating seems 
more pronounced in the higher spatial resolution Model 
B simulation. 

Magnetic field maps are not a good guide to locations 
of large heating; there is much more structure in maps 
showing regions of high heating or maps of heating life¬ 
times than in maps of the magnetic field. Localized heat¬ 
ing events have lifetimes that are about 2 minutes for 
Model A and about 1 minute for Model B. The heating 
events seem to show some recurrence. The intrinsic co¬ 
moving lifetime of events is estimated to be smaller than 
5 minutes. 

In a given location, as a function of time, heating is 
highly episodic and the velocity, density, and tempera¬ 
ture of the plasma will respond on the same timescales as 
well as on the slower timescales dictated by the processes 
cooling the plasma: e.g. enthalpy, flows, conduction. 


and radiative losses. The next step is to construct syn¬ 
thesized observations that clearly link heating events to 
the underlying processes occurring in the chromosphere 
and corona. This may help to clear up the the contra- 
dictory observational results shown by several authors 
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|2002p. A follow up study to pursue 
tween episodic coronal heating and its observational sig¬ 
natures in upper chromospheric, transition region, and 
coronal diagnostics is under way. 


This research was supported by the European Com¬ 
mission funded Research Training Network SOL AIRE. 
It was also supported by the Research Council of Nor¬ 
way through grants of computing time from the Nor¬ 
wegian Programme for Supercomputing and by the 
Swiss National Science Eoundation under grant n*^ 
153302. The research leading to these results has re¬ 
ceived funding from the European Research Council 
under the European Union’s Seventh Eramework Pro¬ 
gramme (EP7/2007-2013) / ERC Grant agreement n^ 
291058 and support by NASA grants NNX08AH45G, 
NNX08BA99G, NNX11AN98G, NNM07AA01C (Hin- 
ode), and NNG09EA40C (IRIS). To analyze the data we 
have used IDL and Vapor (http://www.vapor.ucar.edu). 


REFERENCES 


Abbett, W. P. 2007, ApJ, 665, 1469 

Aschwanden, M. J., Nightingale, R. W., Tarbell, T. D., Wolfson, 

C. J. 2000, ApJ, 535, 1027 

Asgari-Targhi, M., van Ballegooijen, A. A. 2012, ApJ, 746, 81 
Athay, R. G., White, O. R. 1978, ApJ, 226, 1135 
Athay, R. G., White, O. R. 1979, ApJ, 229, 1147 
Baumann, G., Galsgaard, K., Nordlund, A. 2013, Sol.Phys, 284, 
467 

Baumann, G., Haugbplle, T., Nordlund, A. 2012, ArXiv e-prints 
Benz, A. O., Krucker, S. 2002, ApJ, 568, 413 
Berger, M. A., Asgari-Targhi, M. 2009, ApJ, 705, 347 
Berghmans, D., Clette, F., Moses, D. 1998, A& A, 336, 1039 
Bingert, S., Peter, H. 2011, A&A, 530, A112 
Bondeson, A. 1985, Phys. Fluids, 28, 2406 
Bradshaw, S. J., Cargill, P. J. 2010, ApJl, 710, L39 
Cargill, P. J. 1994, ApJ, 422, 381 

Carlsson, M., Hansteen, V. H., de Pontieu, B., McIntosh, S., 
Tarbell, T. D., Shine, D., Tsuneta, S., Katsukawa, Y., Ichimoto, 
K., Suematsu, Y., Shimizu, T., Nagata, S. 2007, PASJ, 59, 663 
Carlsson, M., Leenaarts, J. 2012, A&A, 539, A39 
Carlsson, M., Stein, R. F. 2002, ApJ, 572, 626 
Cavallini, F. 2006, Sol. Phys., 236, 415 
Crosby, N. B., Aschwanden, M. J., Dennis, B. R. 1993, 

Sol. Phys., 143, 275 

Dahlburg, R. B., Klimchuk, J. A., Antiochos, S. K. 2005, ApJ, 
622, 1191 

Datlowe, D. W., Elcan, M. J., Hudson, H. S. 1974, Sol Phys, 39, 
155 

De Pontieu, B., McIntosh, S. W., Carlsson, M., Hansteen, V. H., 
Tarbell, T. D., Boerner, P., Martinez-Sykora, J., Schrijver, 

C. J., Title, A. M. 2011, Science, 331, 55 
De Pontieu, B., McIntosh, S. W., Hansteen, V. H., Schrijver, 

C. J. 2009, ApJL, 701, LI 


De Pontieu, B., Title, A. M., Lemen, J. R., Kushner, G. D., Akin, 

D. J., Allard, B., Berger, T., Boerner, P., Cheung, M., Chou, 
C., Drake, J. F., Duncan, D. W., Freeland, S., Heyman, G. F., 
Hoffman, C., Hurlburt, N. E., Lindgren, R. W., Mathur, D., 
Rehse, R., Sabolish, D., Seguin, R., Schrijver, C. J., Tarbell, 

T. D., Wiilser, J.-P., Wolfson, C. J., Yanari, C., Mudge, J., 
Nguyen-Phuc, N., Timmons, R., van Bezooijen, R., Weingrod, 
L, Brookner, R., Butcher, G., Dougherty, B., Eder, J., 
Knagenhjelm, V., Larsen, S., Mansir, D., Phan, L., Boyle, P., 
Cheimets, P. N., DeLuca, E. E., Golub, L., Gates, R., Hertz, 

E. , McKillop, S., Park, S., Perry, T., Podgorski, W. A., Reeves, 
K., Saar, S., Testa, P., Tian, H., Weber, M., Dunn, C., Eccles, 
S., Jaeggli, S. A., Kankelborg, C. C., Mashburn, K., Pust, N., 
Springer, L., Carvalho, R., Kleint, L., Marmie, J., Mazmanian, 
E., Pereira, T. M. D., Sawyer, S., Strong, J., Worden, S. P., 
Carlsson, M., Hansteen, V. H., Leenaarts, J., Wiesmann, M., 
Aloise, J., Chu, K.-C., Bush, R. L, Scherrer, P. H., Brekke, R, 
Martinez-Sykora, J., Lites, B. W., McIntosh, S. W., 

Uitenbroek, H., Okamoto, T. J., Cummin, M. A., Anker, G., 
Jerram, P., Pool, R, Waltham, N. 2014, Sol. Phys., 289, 2733 

Dennis, B. R. 1985, Sol Phys, 100, 465 
Domingo, V., Fleck, B., Poland, A. 1. 1995, Solphys, 162, 1 
Fossum, A., Carlsson, M. 2005, Nat, 435, 919 
Galsgaard, K., Nordlund, A. 1996, JGR, 101, 13445 
Gudiksen, B. V., Carlsson, M., Hansteen, V. H., Hayek, W., 
Leenaarts, J., Martinez-Sykora, J. 2011, A & A, 531, A154 
Gudiksen, B. V., Nordlund, A. 2005, ApJ, 618, 1020 
Guerreiro, N., Hansteen, V., De Pontieu, B. 2013, ApJ, 769, 47 
Handy, B. N., Acton, L. W., Kankelborg, C. C., Wolfson, C. J., 
Akin, D. J., Bruner, M. E., Caravalho, R., Catura, R. C., 
Chevalier, R., Duncan, D. W., Edwards, C. G., Feinstein, 

C. N., Freeland, S. L., Friedlaender, F. M., Hoffmann, C. H., 
Hurlburt, N. E., Jurcevich, B. K., Katz, N. L., Kelly, G. A., 
Lemen, J. R., Levay, M., Lindgren, R. W., Mathur, D. P., 
Meyer, S. B., Morrison, S. J., Morrison, M. D., Nightingale, 

R. W., Pope, T. P., Rehse, R. A., Schrijver, C. J., Shine, R. A., 
Shing, L., Strong, K. T., Tarbell, T. D., Title, A. M., 

Torgerson, D. D., Golub, L., Bookbinder, J. A., Caldwell, D., 
Cheimets, P. N., Davis, W. N., Deluca, E. E., McMullen, R. A., 
Warren, H. P., Amato, D., Fisher, R., Maldonado, H., 
Parkinson, C. 1999, Solphys, 187, 229 
















20 


Guerreiro et al. 


Hansteen, V. H., Kara, H., De Pontieu, B., Carlsson, M. 2010, 
ApJ, 718, 1070 

Hayek, W., Asplund, M., Carlsson, M., Trampedach, R., Collet, 

R., Cudiksen, B. V., Hansteen, V. H., Leenaarts, J. 2010, A&: 
A, 517, A49 

Hudson, H. S. 1991, SoPhys, 133, 357 
Klimchuk, J. A. 2006, Solphys, 234, 41 
Klimchuk, J. A., Cargill, P. J. 2001, ApJ, 553, 440 
Kosugi, T., Matsuzaki, K., Sakao, T., Shimizu, T., Sone, Y., 
Tachikawa, S., Hashimoto, T., Minesugi, K., Ohnishi, A., 
Yamada, T., Tsuneta, S., Hara, H., Ichimoto, K., Suematsu, 

Y., Shimojo, M., Watanabe, T., Shimada, S., Davis, J. M., Hill, 
L. D., Owens, J. K., Title, A. M., Culhane, J. L., Harra, L. K., 
Doschek, G. A., Golub, L. 2007, Solar Physics, 243, 3 
Krucker, S., Benz, A. O. 1998, ApJL, 501, L213 
Leenaarts, J., Carlsson, M., Hansteen, V., Cudiksen, B. V. 2011, 
A& A, 530, A124 

Lemen, J. R., Title, A. M., Akin, D. J., Boerner, P. F., Chou, C., 
Drake, J. F., Duncan, D. W., Edwards, C. G., Friedlaender, 

F. M., Heyman, G. F., Hurlburt, N. E., Katz, N. L., Kushner, 

G. D., Levay, M., Lindgren, R. W., Mathur, D. P., McFeaters, 
E. L., Mitchell, S., Rehse, R. A., Schrijver, C. J., Springer, 

L. A., Stern, R. A., Tarbell, T. D., Wuelser, J.-P., Wolfson, 

C. J., Yanari, C., Bookbinder, J. A., Cheimets, P. N., Caldwell, 

D. , Deluca, E. E., Gates, R., Golub, L., Park, S., Podgorski, 

W. A., Bush, R. I., Scherrer, P. H., Gummin, M. A., Smith, P., 
Anker, G., Jerram, P., Pool, P., Soufli, R., Windt, D. L., 
Beardsley, S., Clapp, M., Lang, J., Waltham, N. 2012, 

Sol. Phys., 275, 17 
Levine, R. H. 1974, ApJ, 190, 447 

Lin, R. P., Schwartz, R. A., Kane, S. R., Pelling, R. M., Hurley, 
K. C. 1984, ApJ, 283, 421 
Mein, N., Schmieder, B. 1981, A&A, 97, 310 
Nordlund, A. 1982, A&A, 107, 1 

Nordlund, A., Galsgaard, K. 1997, in G. M. Simnett, C. E. 
Alissandrakis,, L. Vlahos (eds.), European Meeting on Solar 
Physics, Vol. 489 of Lecture Notes in Physics, Berlin Springer 
Verlag, 179 

Olluri, K., Gudiksen, B. V., Hansteen, V. H., De Pontieu, B. 

2015, ApJ, 802, 5 
Parker, E. N. 1983, ApJ, 264, 642 
Parker, E. N. 1988, ApJ, 330, 474 

Parnell, C. E., De Moortel, 1. 2012, Royal Society of London 
Philosophical Transactions Series A, 370, 3217 


Parnell, C. E., Jupp, P. E. 2000, ApJ, 529, 554 
Peter, H., Gudiksen, B. V., Nordlund, A. 2006, ApJ, 638, 1086 
Peter, H., Judge, P. G. 1999, ApJ, 522, 1148 
Pontin, D. 1. 2011, Advances in Space Research, 47, 1508 
Priest, E., Forbes, T. 2000, Magnetic reconnection : MHD theory 
and applications / Eric Priest, Terry Forbes. 

Priest, E. R., Demoulin, P. 1995, J. Geophys. Res., 100, 23443 
Priest, E. R., Heyvaerts, J. F., Title, A. M. 2002, ApJ, 576, 533 
Scharmer, G. B. 2006, A&A, 447, 1111 

Scharmer, G. B., Bjelksjo, K., Korhonen, T. K., Lindberg, B., 
Petterson, B. 2003, in S. L. Keil, S. V. Avakyan (eds.). 
Innovative Telescopes and Instrumentation for Solar 
Astrophysics, Vol. 4853 of Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, p. 341 
Scharmer, G. B., Narayan, G., Hillberg, T., de la Cruz Rodriguez, 
J., Lofdahl, M. G., Kiselman, D., Siitterlin, P., van Noort, M., 
Lagg, A. 2008, ApJ, 689, L69 

Schmidt, W., von der Liihe, O., Volkmer, R., Denker, C., Solanki, 

S. K., Balthasar, H., Bello Gonzalez, N., Berkefeld, T., 

Collados, M., Fischer, A., Halbgewachs, C., Heidecke, F., 
Hofmann, A., Kneer, F., Lagg, A., Nicklas, H., Popow, E., 
Puschmann, K. G., Schmidt, D., Sigwarth, M., Sobotka, M., 
Soltau, D., Staude, J., Strassmeier, K. G., Waldmann , T. A. 
2012, Astronomische Nachrichten, 333, 796 

Schrijver, C. J., DeRosa, M. L., Title, A. M. 2010, ApJ, 719, 1083 
Schrijver, C. J., Title, A. M., Berger, T. E., Fletcher, L., 

Hurlburt, N. E., Nightingale, R. W., Shine, R. A., Tarbell, 

T. D., Wolfson, J., Golub, L., Bookbinder, J. A., Deluca, E. E., 
McMullen, R. A., Warren, H. P., Kankelborg, C. C., Handy, 

B. N., de Pontieu, B. 1999, Sol. Phys., 187, 261 
Skartlien, R. 2000, ApJ, 536, 465 
Steffen, M., Muchmore, D. 1988, A&A, 193, 281 
van Ballegooijen, A. A., Asgari-Targhi, M., Cranmer, S. R., 
DeLuca, E. E. 2011, ApJ, 736, 3 
Walsh, R. W., Ireland, J. 2003, A&AR, 12, 1 
Warren, H. P., Brooks, D. H., Winebarger, A. R. 2011, ApJ, 734, 
90 

Welsch, B. T. 2015, PASJ 

Withbroe, G. L., Noyes, R. W. 1977, ArAa, 15, 363 



